Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 25 April 2008 



(MN I^TeX style file v2.2) 



Spatial and Temporal Variations in Interstellar Absorption 
toward HD 72127AB^ 



00 

o 
o 

a., 
< 

(N 

6 



> 

cn 
o 

oo 
O 



Daniel E. Welty^, Thuso Simon^, and L. M. Hobbs^ 

^ University of Chicago, Astronomy and Astrophysics Center, 5640 S. Ellis Ave., Chicago, IL 60637, USA; welty@oddjob.uchicago.edu 
'^University of New Mexico, Department of Physics and Astronomy, Albuquerque, NM 87131, USA; thuso@unm.edu 
^ University of Chicago, Yerkes Observatory, Williams Bay, WI 53191, USA; hobbs@yerkes.uchicago.edu 



25 April 2008 



ABSTRACT 

New optical spectra of Can and Nai toward HD 72127AB provide additional evi- 
dence for both spatial and temporal variations in the complex interstellar absorption 
along the two sight lines; archival UV spectra yield information on the abundances, 
depletions, and physical conditions in the gas toward HD 72127A. Similarities in the 
strengths of various tracers of interstellar material in the two lines of sight suggest 
that the total hydrogen column densities {N ~ 2.5x10^° cm~^) and the depletions 
and ionization in the main components at low LSR velocities also are similar. To- 
ward HD 72127A, the main components are relatively cool (T < 900 K), but with 
depletions resembling those found in warm, diffuse disc clouds; the generally weaker 
components at higher velocities have much milder depletions, more like those found in 
halo clouds. Several trace neutral species - Cai, Cri, and Fei - are much stronger to- 
ward HD 72127B, however. The column density of Cri, for example, is about 30 times 
the value determined toward ^ Oph (the only previous detection of that species in the 
ISM). Dielectronic recombination in warmer gas (T > 5000 K) may be largely respon- 
sible for the enhanced abundances of those trace neutral species toward HD 72127B. 
If the main components toward HD 72127AB are associated with material in the Vela 
SNR, the differences in abundances and physical conditions occur on scales of about 
1100 AU. 

Key words: galaxies: ISM - ISM: abundances - ISM: lines and bands - ISM: struc- 
ture - stars: individual (HD 72127). 
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1 INTRODUCTION 

Small-scale (sub-parsec) spatial structure in the Galactic in- 
terstellar medium (ISM) has been revealed via a variety of 
methods: multi-epoch Hi observations of pulsars (Frail et 
al. 1994; Stanimirovic et al. 2003) , VLBI observations of ex- 
tended extragalactic radio sources (Dieter, Welch, & Rom- 
ney 1976; Brogan et al. 2005), optical spectra of interstellar 
absorption lines toward binary /multiple star systems and 
star clusters (Kemp, Bates, & Lyons 1993; Watson & Meyer 
1996; Points, Lauroesch, & Meyer 2004), and monitoring of 
the optical/UV absorption-line profiles toward stars either 
with significant proper motions or located behind disturbed 



* Based on observations collected at the European Southern Ob- 
servatory, Chile, under program 72.C-0682, and on observations 
made with the NASA/ESA Hubble Space Telescope, which is op- 
erated by the Association of Universities for Research in Astron- 
omy, Inc., under NASA contract NAS5-26555. 
f Visiting observer, European Southern Observatory 



regions (Hobbs et al. 1991; RoUinde et al. 2003; Welty 2007). 
Understanding the nature and origin(s) of that small-scale 
structure in predominantly neutral gas, however, has been 
difficult. Observations of variations in Nal absorption, for 
example, seemed to imply local hydrogen densities (nn) of 
thousands per cm'^ - similar to the values inferred more di- 
rectly from the small-scale variations in H I, but much higher 
than would be expected for clouds in thermal pressure equi- 
librium at typical values uhT 2500 cm~^K (Jenkins & 
Tripp 2001). In principle, the observed small-scale structures 
might be due to geometrical (Heiles 1997), statistical (Desh- 
pande 2000), and/or physical (Lauroesch & Meyer 2003) 
properties of the interstellar clouds. Recent studies combin- 
ing optical and UV absorption-line data have suggested that 
at least some of the variations reflect differences in physical 
conditions (e.g., density, ionization) rather than overall hy- 
drogen column density (Lauroesch et al. 1998; Welty 2007) 
- but that would not explain the variations seen directly 
in A'^(Hl) in other cases. Summaries of recent observational 
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and theoretical studies of small-scale interstellar structure 
may be found in Haverkorn & Goss (2007). 

Observations of stars behind the Vela supernova rem- 
nant (SNR) provided some of the first (and most strik- 
ing) indications of small-scale spatial and temporal variabil- 
ity in interstellar absorption lines (Thackeray 1974; Hobbs, 
Wallerstein, & Hu 1982; Banks & Sembach 1995; Cha & 
Sembach 2000). Such variations - seen toward at least seven 
stars in that region - are not entirely unexpected, given 
the presence of gas at relatively high velocities within the 
SNR (Wallerstein, Silk, & Jenkins 1980; Banks & Sembach 
1995; Cha & Sembach 2000). The spatial and temporal dif- 
ferences in the optical lines of Nal and Call toward the 
binary system HB 72127AB have been particularly notable 
(Thackeray 1974; Hobbs et al. 1982, 1991; Welty, Morton, 
& Hobbs 1996; Cha & Sembach 2000). That system, con- 
sisting of a B2 HI primary [V = 5.20; E{B - V) = 0.10] 
and a B2.5 V secondary [V = 7.09; E{B - V) = 0.10] at 
{l,b) — (262.57,-3.36), is located behind the edge of a fila- 
mentary structure seen in Ha emission, less than a degree 
away from the Vela pulsar. At a distance of about 480 pc, 
however, HB 72127AB is well behind the SNR, which is 
at about 250 pc (Cha, Sembach, & Banks 1999). The 4.5 
arcsec separation between the two stars corresponds to a 
transverse distance of about 2200 AU at the stellar distance 
and about 1100 AU at the distance to the SNR. Thackeray 
(1974) first noted differences in the very broad, very strong 
(>500 mA) interstellar Call K lines toward the two stars; 
Hobbs et al. (1981) then reported both corresponding dif- 
ferences in the Nal B lines and temporal changes in the 
profile of the Can line toward HB 72127A. Betailed fits to 
a series of Can and Nai profiles obtained between 1981 and 
1988 toward HB 72127A indicated continuing variations in 
the column densities and/or velocities of many of the com- 
ponents - particularly those at low LSR velocities (Hobbs 
et al. 1991); high Call/Nal ratios for the higher velocity 
components are suggestive of significant grain disruption. 
Moderate-resolution UV spectra of the two stars obtained 
with the Hubble Space Telescope Goddard High Resolution 
Spectrograph {HST GHRS) revealed both similarities and 
differences in the equivalent widths of the lines from several 
dominant ions; analysis of the C I fine-structure excitation 
implied that some of the gas in both sight lines was charac- 
terized by fairly high thermal pressures (Wallerstein et al. 
1995a; Wallerstein, Vanture, & Jenkins 1995b). 

In this paper, we discuss new optical spectra of both 
HB 72127A and HB 72127B, obtained in 2003 with the ESQ 
VLT Ultraviolet and Visual Echelle Spectrograph (UVES), 
and previously unpublished UV spectra of HB 72127A, ob- 
tained in 1992 with the HST GHRS. The optical and UV 
spectra reveal continuing temporal changes in the interstel- 
lar absorption toward those two stars and provide new in- 
formation on the abundances and physical conditions in the 
gas. In Section [2] we describe the observed spectra and our 
reduction and analysis techniques. In Section |3] we discuss 
the elemental abundances and depletions in the gas, the 
temporal changes in various features, the differences in ab- 
sorption between the two sight lines, and the remarkably 
strong absorption from some trace neutral species seen to- 
ward HB 72127B. In Section |4l we provide a summary of 
our results. 



2 DATA 

2.1 Observations and data processing 

Optical spectra of HB 72127A and HB 72127B were ob- 
tained on 2003 November 14 using the ESO/VLT UT2 tele- 
scope and UVES spectrograph (Bekker et al. 2000), during 
a run devoted primarily to observations of weak atomic and 
molecular interstellar lines in the Magellanic Clouds (Welty 
et al. 2006; Welty & Crowther, in preparation). The stan- 
dard dichroic 1 390/564 setting and a slit width correspond- 
ing to 0.7 arcsec were used to obtain nearly complete cov- 
erage of the wavelength range 3260-6680 A (on three CCB 
detectors: two EEV CCB-44 and one MIT/LL CCIB-20) at 
resolutions of about 4.5 km s~^ in the blue and about 4.9 
km s^^ in the red. This setup includes lines from Nal (U 
and B doublets at 3302 and 5889/5895 A), Cai (4226 A), 
Can (3933/3968 A), and Tin (3383 A); the strongest lines 
from CH (4300 A), CH+ (4232 A), and ON (3874 A); and a 
number of the diffuse interstellar bands. Standard routines 
within IRAF were used to remove bias and to divide sec- 
tions of the 2-B images containing the spectral order(s) of 
interest by a normalized fiat-field derived from quartz lamp 
exposures. The 1-B spectra then were extracted from the 
fiat-fielded image segments via the APEXTRACT routines, us- 
ing variance weighting (with the appropriate values for read 
noise and gain for each detector). Wavelength calibration 
was accomplished via Th-Ar lamp exposures, which were 
obtained at the beginning and end of each night, using the 
thorium rest wavelengths tabulated by Palmer & Engelman 
(1983). The spectra were then normalized via Legendre poly- 
nomial fits to the continuum regions surrounding the inter- 
stellar (and stellar) absorption lines. Equivalent widths for 
the various interstellar absorption features were measured 
from those normalized spectra. The empirical signal-to-noise 
(S/N) ratios measured in the continuum regions are typi- 
cally 100-200 per half resolution element - yielding 3a 
uncertainties or detection limits for weak, unresolved ab- 
sorption lines of 1-2 mA. Equivalent widths measured for 
various interstellar lines toward several other Galactic stars 
observed during the run show very good agreement with pre- 
viously reported values (Gredel et al. 1991, 1993; Albert et 
al. 1993; Welty et al. 1999b). Observations of the bright, 
nearby, rapidly rotating star tp'^ Aqr were used to check 
for weak atmospheric absorption features and instrumental 
artefacts near the interstellar lines of interest. 

UV spectra of HB 72127A were obtained in 1992 April 
and November with the HST GHRS, under GO programs 
2251 and 3993 (L. M. Hobbs, PI). The data set includes 
13 high-resolution (FWHM = 3.3-3.7 km s"^) ECH-B set- 
tings, with central wavelengths from 1745 to 2800 A, and ten 
medium-resolution (FWHM = 13-20 km s'^) G160M set- 
tings, with central wavelengths from 1135 to 1664 A. Multi- 
ple, slightly offset sub-exposures were obtained at each set- 
ting - using the WSCAN procedure for the ECH-B spectra 
and the FP-SPLIT procedure for the G160M spectra - in 
order to identify and reduce the effects of detector fixed- 
pattern noise. The individual sub-exposures for each wave- 
length setting were aligned via fits to the interstellar absorp- 
tion features, cleaned, and co-added; the summed spectra 
were then normalized via Legendre polynomial fits to the 
continuum regions. The empirical S/N ratios in the normal- 
ized spectra range from about 35 to 100 per half resolution 
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Figure 1. Interstellar Nal and Call profiles toward HD 72127A 
and HD 72127B. For each profile (or pair of profiles), the sight 
line, species, observation date, and resolution are given at the 
right. The UVES spectra (2003 Nov) reveal striking differences 
between the profiles for the two lines of sight, which are separated 
by only 4.5 arcsec. Several earlier profiles toward HD 72127A 
(Hobbs et al. 1991; Welty et al. 1996) are included to show tem- 
poral variations. Most of the components in Call exhibit some 
change over the 15-year span represented here. 



element for the ECH-B data and from about 80 to 150 for 
the lower resolution G160M data. More detailed discussions 
of the GHRS data and our adopted reduction and analysis 
procedures may be found in Welty et al. (1999b). 

2.2 Spectra and equivalent widths 

The normalized line profiles for Call and Nal toward 
HD 72127A and HD 72127B are shown in Fig. [T] In each 
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Figure 2. Profiles of weaker interstellar lines toward HD 72127A 
(upper for each pair) and HD 72127B (lower for each pair) ob- 
served with UVES. The weaker member of the Nal A3302 doublet 
is seen near -|-70 km s~^. The lines from Nal and Tin are similar 
in strength toward the two stars, but the lines from Ca I, Cr I, and 
Fel are much stronger toward HD 72127B. 



case, the stronger member of the doublet (if present) is 
shown at a continuum level of 1.0, and the weaker member 
of the doublet (if present) is offset by -1-0.2. Velocities with 
respect to the local standard of rest (LSR) may be obtained 
by subtracting 13.4 km s~^ from the heliocentric velocities 
shown in the figure^. Several previously observed profiles for 
HD 72127A - Can and Nai from 1988 (FWHM ~ 2.5-2.7 
km s"^; Hobbs et al. 1991) and Call from 1994 (FWHM ~ 
0.3 km s~^ Welty et al. 1996) - are included for compari- 
son. Spectra of some of the weaker optical lines are shown 
in Fig. m where the profile for HD 72127A is offset above 
the one for HD 72127B in each case (and note the expanded 
vertical scale for all the lines) . Profiles of some of the U V ab- 
sorption lines toward HD 72127A are shown in Fig.|3l where 
(again) the vertical scale has been expanded for the weaker 
lines. Equivalent widths for the various lines measured from 
the normalized optical and UV spectra are listed in Table [1] 
for comparison, values reported for Mgll, Sll, and Cell by 
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Table 1. Equivalent widths. 
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Figure 3. Selected UV absorption lines toward HD 72127A, ob- 
tained in 1992 with the HST GHRS, at resolutions of about 3.3- 
3.7 km s^-*^. The absorption near 60 km s^^ in the Znll spectrum 
is due to Mgi. Note the expanded vertical scale for several of the 
weciker lines. 



Species 


A 




TJT\ vol OV A 




Oi 


1355.5977 


-2.803 


<2.6 




Nai 


3302.3690 


1.472 


2.1±0.9 


■ 




3302.9780 


1.170 


1.9±0.8 


2.7±1.0 




5889.9510 


3.586 


216.0±4.2 


277.4±4.2 




5895.9242 


3.285 


135.9±4.0 


187.6±3.3 


Mgi 


1827.9351 


1.646 


4.4±0.8 






2026.4768 


2.360 


18.3±1.3 




Mgii 


1239.9253 


-0.106 


14.8±1.7 










[14.9] 


[13.0] 




1240.3947 


-0.355 


9.5±1.2 










[8.6] 


[8.0] 

L J 


All 


3944.0060 


2.664 


<1.1 


<l.l 


Aim 


1854.7184 


3.016 


169.1±2.6 






1862.7910 


2.714 


100.6±2.2 




Si I 


2515.0725 


2.725 


1.8±0.4 




Si II 


1808.0129 


0.575 


115.1±4.1 




So 


1250.5780 


0.832 


[111.5] 


[104.0] 




1253.8050 


1.136 


[177.1] 


[177.1] 




1259.5180 


1.320 


[209.0] 


[204.3] 


Oil 


1347.2396 


2.314 


6.2±2.2 




Cai 


4226.7280 


3.874 


1.0±0.4 


16.0±0.9 


Can 


3933.6614 


3.392 


531.8±1.9 


508.4±2.4 




3968.4673 


3.092 


327.4±2.0 


339.2±3.6 


Til 


3635.4620 


2.962 


<2.0 


<1.6 


Tin 


3383.7588 


3.083 


17.0±2.6 


21.3±4.2 


On 


3578.6840 


3.117 


<1.7 


5.4±0.7 




3593.4810 


3.019 




5.5±0.6 




3605.3210 


2.911 




4.5±0.6 




4254.3320 


2.670 




2.7±0.4 




4274.7960 


2.556 




1.8±0.4 




4289.7160 


2.427 




1.2±0.4 


Crii 


2056.2569 


2.326 


22.4±4.5 






2062.2361 


2.195 


17.7±4.0 




Mm 


4030.7530 


2.357 


<1.2 


<1.2 


Fei 


2484.0209 


3.131 


1.7±0.5 






2523.6083 


2.710 


0.9±0.3 






3440.6057 


1.910 




1.3±0.6 




3719.9347 


2.184 


<1.3 


3.5±0.7 




3859.9114 


1.923 




2.1±0.4 


Fell 


2249.8768 


0.612 


35.0±3.1 






2374.4612 


1.871 


384.0±3.3 




Nil 


3369.5656 


1.913 


<2.2 


<2.6 


Nili 


1741.5531 


1.871 


32.3±5.5 




Znii 


2026.1370 


3.007 


56.3±3.6 






2062.6604 


2.705 


42.6±2.9 




Go II 


1237.0589 


3.183 


2.9±1.1 










4-(l 





Entries are equivalent width ± la uncertainty (niA); limits are 
3iT. Wavelengths (in air above 3000 A) and /-values are from 
Morton (2003). Values for lines above 3000 A are from 2003 
UVES spectra; values for lines below 3000 A are from 1992 
GHRS spectra; values in square brackets are from 1992-1993 
GHRS spectra (Wallerstein et al. 1995). 



Wallerstein ct al. (1995b) (from GHRS spectra obtained in 
1992 and 1993) arc given in square brackets. 

In both sight lines, the trace neutral species (e.g., Nai, 



^ The conversion to LSR velocities assumes a solar motion of 20 
km s-i toward {a,S) = (181^ ,30°). 
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Mgl) are found primarily in an apparently single compo- 
nent at +10 to +13 km (ulsr ~ ^3 to km s~^), 
with much weaker absorption extending from about —20 to 
+55 km s~^. Toward HD 72127A, the absorption from dom- 
inant ions of typically mildly depleted elements (e.g., Znll) 
is concentrated in several components between about +10 
and +20 km s~^, with hints of much weaker absorption at 
lower and higher velocities. The absorption from dominant 
ions of typically more severely depleted elements (e.g., Sin, 
Fell) also is strongest between +10 and +20 km s~^, but 
becomes increasingly prominent at both lower and higher 
velocities for the species most affected by depletion. Can 
- which can be affected by both ionization and depletion 
(e.g., Welty et al. 1996) - shows moderately strong absorp- 
tion over the entire velocity range. 

2.3 Profile fits 

Multicomponent fits to observed complex absorption-line 
profiles enable blended lines and components to be disentan- 
gled, account explicitly for saturation effects in the stronger 
components, and allow spectra obtained at (somewhat) dif- 
ferent resolutions to be compared. The method of profile 
fitting was therefore used to estimate column densities (A*'), 
line widths (fe ~ FWHM/1.665), and velocities {v) for the in- 
dividual components discernible in the spectra (e.g., Hobbs 
et al. 1991; Welty, Hobbs, & Morton 2003). For each species, 
the minimum number of components needed to achieve 'sat- 
isfactory' fits to the detected lines was adopted (given the 
resolution and S/N ratios characterizing the spectra and as- 
suming that each component may be represented by a sym- 
metric Voigt profile). Multiple lines for a given species (e.g., 
the two lines of Ca n) were fitted simultaneously. The struc- 
ture determined for Can was used to estimate individual 
component column densities for Tin from the weak A3383 
lines. The lines from Sin, Crn, Fen, Nin, and Znn seen in 
the GHRS ECH-B spectra were fitted simultaneously - using 
a common velocity structure but adjusting the component 
column densities (and some &- values) for each species. 

The adopted interstellar component parameters derived 
from the UVES spectra of HD 72127A and HD 72127B 
are listed in Table [J] the component parameters obtained 
from the simultaneous fit to the GHRS ECH-B spectra of 
HD 72127A are given in Table |3l For each component in 
each sight line. Table [2] gives the component number, the 
heliocentric velocity and &- value (for Nal and Can), and the 
column densities of Nal, Cal, Can, and Tin. Component 
column densities or 6- values in square brackets were fixed, 
but relatively well determined in the fits; values in paren- 
theses also were fixed, but are less well determined. For Nal, 
the column density of the strongest component in each sight 
line was fixed at a value consistent with that derived from 
the weak A3302 doublet; the parameters for the other com- 
ponents were derived from fits to the much stronger Nal D 
lines. Independent fits to the profiles of the Na I D lines and 
the Can H and K lines yielded rather similar component 
velocities (in both sight lines), and components seen in the 
two species are considered to be associated if the velocities 
agree within about 1-2 km s~^. 

The 1 G uncertainties given for A'', 6, and v in Tables [2] 
and [3] reflect only photon noise, and thus should be consid- 
ered as lower limits to the true uncertainties. Consideration 



Table 4. Column densities. 



Species HD 72127A HD 72127B 23 Ori SLV C Oph B 

Htot [20.40] [20.40] 20.71 21.15 

Oi <17.14 - 17.29 17.63 

Nai 12.46+0.12 12.66+0.06 13.41 13.85 

Mgi 12.92+0.04 - 13.80 13.80 

Mgii 15.29+0.04 - 15.57 15.32 

Sii 11.20+0.11 - 11.79 

Sin 15.43+0.01 - 15.37 15.34 

All <10.83 <10.83 <11.53 <10.22 

CI I 12.52+0.16 - 13.29 14.48 

Cai 9.58+0.08 10.76+0.02 10.20 9.46 

Can 12.99+0.01 13.22+0.02 12.04 11.71 

Til <10.88 <10.79 <10.72 

Tin 11.68+0.06 11.78+0.08 11.15 11.03 

Cri <10.61 11.26+0.02 <10.59 9.78 

Crn 12.80+0.04 - 12.64 12.39 

Mm <11.17 <11.17 10.23 <11.30 

Mno - - 13.08 13.20 

Fei 10.77+0.11 11.89+0.03 11.34 11.23 

Fell 14.69+0.01 - 14.30 14.32 

Nil <11.95 <12.03 <10.26 <11.40 

Nil! 13.48+0.03 - 12.90 12.91 

Znn 12.72+0.01 - 13.24 

Gcii 11.25+0.17 - 11.86 11.75 



Entries are log [A'' (cm~^)] + la uncertainty, based on fits to 
the line profiles; limits are 3(T. Values for 23 Ori SLV (main) 
components are from Welty et al. (1999b), except for Nal, Til, 
and Cri, which are from new UVES data. Values for C, Oph B 
(main) component are from Welty &; Hobbs (2001), Welty et al. 
(2003), and references therein. 



of the effects of uncertainties in continuum placement (most 
significant for the stronger lines) and of allowing fixed pa- 
rameters (usually h) to vary would increase the uncertain- 
ties in the various parameters (primarily N). In view of the 
statistics of component width and separation obtained from 
higher resolution, higher S/N spectra of Galactic targets 
(Welty et al. 1994, 1996; Welty & Hobbs 2001), however, it 
is likely that the 'true' component structures are even more 
complex than those listed in the tables - and that uncer- 
tainty may well be even more significant for the individual 
component parameters. Differences in the component veloc- 
ities derived for Nal and Can, for example, may be due to 
component-to-component differences in relative abundances 
in an underlying more complex structure. Higher resolution 
spectra of Can toward HD 72127A do exhibit more complex 
structure (Fig.[T] Welty et al. 1996). Differences between the 
component velocities derived from optical and UV spectra 
of HD 72127A may refiect (in addition) possible temporal 
variations between 1992 and 2003. 

The total sight line column densities are given in Ta- 
bled with values for the main neutral components toward 
23 Ori (Wehy et al. 1999b) and C Oph (Wehy et al. 2003 
and references therein) given for comparison. In all cases, the 
values derived from the profile fits are consistent with esti- 
mates or limits based on integrating the 'apparent' optical 
depth over the absorption lines; the total sight line column 
densities are in most cases fairly well determined. 



[20.40] 




[20.40] 


20.71 


21.15 


<17.14 




- 


17.29 


17.63 


12.46+0.12 


12, 


,66+0.06 


13.41 


13.85 


12.92+0.04 




- 


13.80 


13.80 


15.29+0.04 




- 


15.57 


15.32 


11.20+0.11 




- 


11.79 


- 


15.43+0.01 




- 


15.37 


15.34 


<10.83 




<10.83 


<11.53 


<10.22 


12.52+0.16 




- 


13.29 


14.48 


9.58+0.08 


10, 


,76+0.02 


10.20 


9.46 


12.99+0.01 


13, 


,22+0.02 


12.04 


11.71 


<10.88 




<10.79 


<10.72 




11.68+0.06 


11, 


,78+0.08 


11.15 


11.03 


<10.61 


11, 


,26+0.02 


<10.59 


9.78 


12.80+0.04 






12.64 


12.39 


<11.17 




<11.17 


10.23 


<11.30 








13.08 


13.20 


10.77+0.11 


11, 


,89+0.03 


11.34 


11.23 


14.69+0.01 






14.30 


14.32 


<11.95 




<12.03 


<10.26 


<11.40 


13.48+0.03 






12.90 


12.91 


12.72+0.01 






13.24 




11.25+0.17 






11.86 


11.75 
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Table 2. Component structures (optical data). 







Nal 




Cai 




Ca II 




Till 


Caii/Nai 


Caii/Tiii 


Comp 


V 


b 


Nio 


Ng 


V 


6 


ATio 










(km s-l) 


(km s-i) 


(cm-2) 


(cm-2) 


(km s-l) 


(km s-i) 


(cm-2) 


(cm-2) 






HD 72127A 


1 










-27.7±2.1 


(3.0) 


1.0±0.5 








2 


-13.8±0.3 


4.4±0.5 


9.8±0.6 




-15.4±0.1 


3.2±0.1 


218.3±3.8 


2.6±1.4 


22.3 


~84.0 


3 










-8.2±0.1 


(4.0) 


112.6±3.0 


2.0±1.5 




~56.3 


4 


-l.litO.6 


(3.0) 


8.4±1.2 
















5 


3.7±0.2 


(2.0) 


21.3±1.4 




2.9±0.1 


3.3±0.2 


120.1±3.0 


4.7±1.3 


5.6 


25.6 


6 


12.8±0.1 


1.2±0.1 


[250.0±70.] 


3.8±0.7 


12.6±0.2 


3.6±0.4 


148.4±12. 


15.9±1.6 


0.6 


9.3 


7 


18.8±0.2 


(3.0) 


14.0±0.8 




18.3±0.5 


(3.0) 


41.3±9.4 


13.6±1.5 


3.0 


3.0 


8 










26.0±0.2 


(3.0) 


67.1±1.8 


2.8±1.4 




~24.0 


9 


32.8±0.5 


(4.0) 


7.4±0.7 




32.3±0.2 


(3.0) 


81.6±2.2 


3.3±1.4 


11.0 


~24.7 


10 


40.7±1.3 


(4.0) 


3.0±0.6 




40.0±0.1 


3.8±0.1 


130.2±2.9 


3.0±1.4 


43.4 


~43.4 


11 


52.0±1.3 


(5.0) 


2.5±0.5 




54.7±0.1 


4.9±0.1 


61.6±1.0 




24.6 




HD 72127B 


1 


-15.1±1.3 


(2.0) 


1.0±0.4 




-16.0±0.1 


3.5±0.1 


.58.1±1.1 




58.4 




2 










-5.4±0.1 


2.5±0.2 


.51.7±1.2 


1.7±1.2 




~30.4 


3 


-1.1±0.4 


(2.5) 


6.9±0.7 
















4 


6.5±0.4 


(3.0) 


31.7±4.3 




4.3±0.3 


(3.0) 


113.1±17. 




3.6 




5 


12.2±0.1 


2.5±0.1 


[400.0±55.] 


57.9±1.2 


10.5±0.1 


(4.0) 


1060.0±65. 


28.4±1.9 


2.7 


37.3 


6 


21.4±0.2 


(3.5) 


18.2±0.7 


5.5±0.9 


21.0±0.1 


(4.0) 


185.7±3.0 


20.2±1.5 


10.2 


9.2 


7 










29.4±0.2 


(3.0) 


41.2±1.5 








8 


34.7±0.4 


(3.0) 


5.7±0.5 




35.7±0.3 


(3.0) 


27.2±1.3 




4.8 




9 


43.3±1.2 


(4.0) 


2.7±0.5 




42.9±0.2 


(3.0) 


36.9±1.2 




13.7 




10 


54.2±1.2 


(4.0) 


2.1±0.5 




.53.8±0.1 


5.6±0.1 


90.4±1.6 


3.5±1.5 


43.0 


25.8 


11 










88.0±1.0 


(4.0) 


4.0±0.8 








12 










96.9±1.3 


(4.0) 


3.3±0.8 








The subscript on each column density header indicates the power of ten by which the listed value of that species should be multiplied. 


Toward HD 72127B, Cr i and Fe I are 


detected, with [AT, v] equal to [(1.8±0.1)xl0", ll.OiO.l] and [(7.8±0.6)xl0", 


10.6±0.2], respectively. 


Table 3. Component structures (HD 72127A UV data). 














Comp 


V 


b 


N(Zn II) 


iV(Siii) 


Ar(Crii) 


N{Fe II) 


Ar(Niii) 


Feii/Zn II 






(km s-l) 


(km s-i) 


(lOlO cm-2) 


(10^3 cm- 


2) (IQii cm-2) (^0^2 cm-2) (^^^12 ^^-2^ 












9.0 


3.5 


4.5 




2.5 


2.0 






1 


-24.3±1.2 


(3.0) 










1.1±0.3 








2 


-15.9±0.1 


2.5 




6.6±0.9 6.0±1.3 


59.8±2.8 


2.9±0.6 


>665 




3 


-8.5±0.2 


4.0 




11.8±1.1 6.8±1.5 


52.7±2.4 


4.8±0.7 


>585 




4 


2.9±0.1 


2.8 


17.3±3.8 


24.6±1.3 5.8±1.3 


54.6±1.4 


3.4±0.6 


315 




5 


13.2±0.2 


[3.7] 


338.2±13. 


111.7±5.1 26.7±1.9 


115.1±7.7 


6.5±0.7 


35 




6 


19.5±0.3 


[4.0] 


150.5±7.7 


58.2±3.5 9.2±2.1 


85.1±1.7 


3.9±0.7 


55 




7 


28.0±0.1 


2.4 




7.7±1.1 3.9±1.4 


31.1±1.7 


2.2±0.6 


>345 




8 


32.9±0.5 


(3.5) 




13.7±1.3 




15.1±1.7 




>170 




9 


39.9±0.1 


3.4 


5.0±3.9 


12.0±1 


.8 4.4±1.4 


37.4±1.0 


3.5±0.6 


750 




10 


44.6±1.2 


(3.5) 




5.5±1.4 












11 


54.7±0.2 


5.7 




16.6±2.0 




32.9±1.2 


2.5±0.7 


>365 




12 


57.8±0.8 


(2.0) 










2.3±1.0 










10-20 km s 


-1 


488.7±15. 


169.9±6.2 35.9±2.8 


200.2±7.9 


10.4±1.0 


40 






outliers 




<71.0 


98.1±5.1 <36.6 


286.8±7.7 


19.5±1.7 


>405 






total 




(526.0±18.) 


268.0±8.0 (63.0±5.0) 


487.0±11. 


29.9±2.0 


95 





Total column densities for Znii and Cm in parentheses should be considered as lower limits, as not all components were detected. 
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Table 5. Abundance and depletion data (HD 72127A). 



3.1 Abundances and depletions 

In view of the spectral types of HD 72127A and HD 72127B, 
it is likely that the Lyman-Q absorption observed toward the 
two stars (in archival /(7iJ spectra) is an inseparable blend of 
stellar and interstellar components. Several indirect meth- 
ods thus have been used to estimate the interstellar hydro- 
gen column densities toward the two stars: 

(1) The first makes use of the strong, nearly linear corre- 
lation observed between the equivalent width of the diffuse 
interstellar band (DIB) at 5780 A and iV(Hl) (e.g., Herbig 
1993, 1995; Welty et al. 2006; D. G. York et al., in prepa- 
ration). The equivalent widths estimated from our UVES 
spectra - W(5780) = 15±6 mA for HD 72127A and W(5780) 
= 12±5 mA for HD 72127B - suggest log[7V(H l)] of about 
20.2 for HD 72127A and about 20.1 for HD 72127B. The 
observed W(5780) are somewhat low, relative to the corre- 
sponding A^(Nal), however, suggesting that the 5780 A DIB 
may be weaker than usual toward HD 72127AB - as has been 
observed, for example, toward some Galactic and Magellanic 
Clouds targets (Herbig 1993; Welty et al. 2006). The iV(Hl) 
estimated from W(5780) might thus be best viewed as lower 
limits. 

(2) An estimate for N{B) toward HD 72127A may be ob- 
tained from the observed column density of Znll, as zinc 
is generally only very mildly depleted in diffuse interstellar 
clouds. If the depletion of zinc D(Zn) = —0.2 dex (typical for 
warm diffuse gas), then log[A'^(Znll)] = 12.72 implies that 
log[Af(H)] would be about 20.3, for a solar zinc abundance 
of —7.37 dex (Lodders 2003). Fits to the lower resolution 
G160M spectra of the Sll lines at 1250, 1253, and 1259 A 
using the component structures found for Zn ll and Si ll yield 
consistent estimates for A^(H), assuming a solar abundance 
of —4.81 dex and no depletion for sulphur. 

(3) If the typical Galactic gas-to-dust ratio log[7V(H)/i5(_B — 
V")] = 21.71 (e.g., Welty et al. 2006) applies to the sight line 
toward HD 72127AB, then the observed E{B - V) = 0.10 
would imply log[iV(H)] ~ 20.7±0.3 toward both stars. 

(4) The upper limit on A''(0 1) (from the non-detection of 
the weak A1355 line) toward HD 72127A may be used to 
place a corresponding upper limit on A*' (Hi), as the two 
species are strongly coupled via charge exchange. The ob- 
served log[Ar(Ol)] < 17.14 implies log[A'(Hl)] < 20.55, for 
a solar oxygen abundance of —3.31 dex and a typical 'warm 
cloud' depletion of —0.1 dex (e.g., Cartledge et al. 2004). 

(5) The equivalent widths observed for the weak lines of 
typically moderately depleted Mgll (1239 and 1240 A) and 
for the strong lines of typically undepleted Sll (1250, 1253, 
and 1259 A) are very similar toward the two stars (Table [T] 
Wallerstein et al. 1995b) - suggesting that A''(H) also is sim- 
ilar toward both. 

(6) The upper limit on log[Af(CH)] < 12.4 obtained from 
the UVES spectra suggests that log[A'^(H2)] is less than 20.0 
toward both stars (e.g., Welty et al. 2006). 

As the different indicators (photometry, optical spectra, UV 
spectra) are from different epochs, these comparisons and 
estimates are somewhat uncertain, given the observed tem- 
poral variations for some interstellar quantities in both sight 
lines. None the less, in view of the various estimates and 
limits, we have adopted log[7V(Hl)] ~ log[A'(H)] ~ 20.4, 



Element 


Solar 


D(Warm)/ 
D(Cold) 


main 


D(72127A) 
outlier 


total 


O 


8.69 


-0.1/-0.2 


<0.03 




< -0.05 


Na 


6.30 


-0.6/-0.6 








Mg 


7.55 


-0.6/-1.3 






-0.66 


Si 


7.54 


-0.4/-1.3 


-0.63 


-0.15 


-0.51 


Ca 


6.34 


-2.0/-3.6 


-2.38 




-1.75 


Ti 


4.92 


-1.3/-2.9 


-1.77 


-1.25 


-1.64 


Cr 


5.65 


-1.2/-2.2 


-1.41 


-0.82 


-1.25 


Fe 


7.47 


-1.3/-2.2 


-1.49 


-0.61 


-1.18 


Ni 


6.22 


-1.4/-2.3 


-1.52 


-0.53 


-1.14 


Zn 


4.63 


-0.2/-0.4 


-0.26 




-0.31 


Ge 


3.62 


-0.6/-0.9 


-0.69 







Solar abundances are from Lodders (2003). Representative warm 
and cold cloud depletions are updates of values in Welty et al. 
(1999b). Last three columns give depletions for HD 72127A: 
main components (10-20 km s~^), all other components, and all 
components, respectively. Those depletions assume log[Af(H)] = 
20.32, 19.60, and 20.40 dex for the three groups. 

with an uncertainty of perhaps ±50 percent, toward both 
HD 72127A and HD 72127B. 

The fits to the Znll, Sll, and Sill profiles toward 
HD 72127A suggest that roughly 80-90 percent of the total 
Hi is found in the components between 10 and 20 km s~^ 
('main components'), with the remainder in the outlying 
higher velocity components. On the other hand, roughly 
60 percent of the Fell and Nill is found in the outlying 
components. Inspection of the ratios of the column densi- 
ties of typically depleted species to that of (nearly) unde- 
pleted Znll (e.g., Fell/Znil; last column of Table [3)1 indi- 
cates that the depletions in the outlying components are 
significantly less severe than those in the main components. 
If the adopted total A'^(H) is apportioned to the main com- 
ponents (2.1 X lO'^" cm~^) and the outlying components (4 
X 10^^ cm~^), then the depletions may be estimated sep- 
arately for those two component groups (Table (5)1. For the 
main components, the resulting depletions are broadly con- 
sistent (within the uncertainties) with those in 'warm dif- 
fuse' clouds in the Galactic disc (Savage & Sembach 1996; 
Welty et al. 1999); for the outlying components, the deple- 
tions are more similar to those in Galactic halo clouds. The 
similarities in E{B — V), A'^(Till), and the equivalent widths 
of the 5780 A DIB and the lines of Mg ll and S ll for both 
HD 72127A and HD 72127B (Wallerstein et al. 1995b; this 
paper) suggest that both the total Ai'(H) and the average de- 
pletions are also similar (within factors of about 2) for the 
two lines of sight; the depletions may be slightly less severe 
toward HD 72127B. 

Examination of the profiles of the Aim lines at 1854 
and 1862 A may provide some insight into the distribution 
of ionized gas toward HD 72127A. Integrating the appar- 
ent optical depths over the profiles yields a total Af(Allll) of 
about 1.7x 10^^ cm~'^ in each case, so the lines do not appear 
to be significantly saturated. There is very little Aim ab- 
sorption at the velocities of the main components between 
10 and 20 km s"'^ - which thus are likely to be predom- 
inantly neutral. The Aim absorption at higher and lower 
velocities (especially the strong absorption centred at about 
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35 km s~^) signals the presence of some ionized gas at those 
outlying velocities. Crude fits to the Al ll line at 1670 A (us- 
ing the Fell component structure), however, suggest that the 
Alll/Allll ratio is about 3-4 (on average) for those outlying 
components - so that the fraction of ionized gas may be rel- 
atively small there as well. Uncertainties as to the amounts 
of Hll and Xlll (i.e., the second ions of elements whose first 
ions Xll are dominant in Hi gas) present in the outlying 
components imply corresponding uncertainties in the deple- 
tions estimated for those components in Table [5] 

3.2 Call/Nal ratios 

The Call/Nal ratio is usually considered to be an indicator 
of the (highly variable) depletion of calcium, and can range 
from < 0.01 in cold, dense clouds (where calcium generally is 
severely depleted) to > 10 in warm, diffuse gas (where some 
fraction of the calcium has been returned to the gas phase). 
Hobbs et al. (1982) found unusually high Call/Nal ratios 
(>9) for all five of the components toward HD 72127A seen 
in 1977. The detailed fits reported by Hobbs et al. (1991), 
however, yielded average ratios ranging from 1.4 for the com- 
ponent at -1-13 km s~^ (with strongest Nal) to >30 for the 
components at —19 and —8 km s~^, for observations be- 
tween 1981 and 1988. Fits to the 2003 UVES spectra imply 
Call/Nal ratios ranging from 0.6 (at -|-13 km s"'^) to >20 
(at -15, -f40, and -1-55 km s"^) toward HD 72127A and 
from 2.7 (at -1-11 km s"^) to >40 (at -16 and -1-54 km s'^) 
toward HD 72127B. At least some of the differences between 
the current and former results are likely due to the more de- 
tailed and precise component information obtained from the 
more recent spectra and to the availability of A3302 data 
for Nal. The depletions that would be inferred from the 
Call/Nal ratios in the different components thus are con- 
sistent with the depletion levels determined more directly 
from the various dominant species in the previous section. 

3.3 Temporal variations 

Temporal variations over the past 30 years in the complex 
profiles of Nal and Call toward HD 72127A have been well 
documented. Hobbs et al. (1982) noted the appearance of 
a 'new' component at -1-15 km s~^ in a Call spectrum ob- 
tained in 1981 April. Detailed fits to that same spectrum 
(and others obtained through 1988 Dec) revealed both ad- 
ditional weaker components and variations in the column 
densities and/or velocities of most of the Call components 
between —15 and -f 17 km s~^ (Hobbs et al. 1991). An ultra- 
high resolution Ca ll spectrum obtained in 1994 Feb revealed 
continued variations over that velocity range, more complex 
structure between —20 and -1-40 km s~^, and new varia- 
tions between -1-25 and -1-35 km s~^ (Fig. [1] Welty et al. 
1996; see also Cha & Sembach 2000). The UVES Can spec- 
trum (2003 Nov) shown in Fig. [T] indicates further changes 
in the component column densities over most of the veloc- 
ity range and continued evolution of the structure between 
-1-25 and -1-35 km s~^. The component velocities found for 
the UVES spectrum are very similar, however, to those de- 
termined from spectra of comparable resolution obtained in 
1994 and 1996 (Cha & Sembach 2000), when allowance is 
made for differences in adopted component structure (Ta- 
ble [6]). The most noticeable change in the Nal profile toward 



Table 6. Can component velocities (1988-2003). 



1988 Dec 


1994 Feb 


1994 Feb 


1996 J/F 


2003 Nov 


HD 72127A 










-27.7 


-18.8 


-17.7 








— 14.6 


-15.3 
-9.1 


-15 


-15 


— 15.4 


—8.0 


-8.2 
-0.8 


-8 


-8 


-8.2 




2.0 


3 


3 


2.9 


3.5 


4.6 
7.8 








11.0 


12.7 
13.2 






12.6 


17.4 


14.3 
22.1 


14 


14 


18.3 


27.7 


29.8 
31.1 
38.5 


30 


29 


26.0 
32.3 


39.8 


40.9 


41 


40 


40.0 


54.8 


53.9 


55 


55 


54.7 


59.8 


58.2 









HD 72127B 



-15 


-16.0 


-6 


-5.4 




4.3 


13 


10.5 




21.0 


30 


29.4 




35.7 


40 


42.9 


55 


53.8 




88.0 




96.9 



Velocities are heliocentric; ultra-high resolution spectrum from 
1994 Feb used as standard for HD 72127A. References: 1988 Dec 
(Hobbs et al. 1991); 1994 Feb (Welty et al. 1996); 1994 Feb and 
1996 J/F (Cha & Sembach 2000); 2003 Nov (this paper). See 
Hobbs et al. (1982, 1991) for values from previous epochs. 



HD 72127A since 1988 is the appearance of a distinct com- 
ponent near +4 km s~^ (Cha & Sembach 2000; this paper) 
superimposed on (or in place of) the smooth, broad shoul- 
der seen in previous (higher resolution) spectra short-ward 
of the strongest Nal absorption near -1-13 km s~^. 

Unfortunately, much less information is available for 
the interstellar lines toward the fainter HD 72127B. Com- 
ponent velocities for Can reported by Hobbs et al. (1982) 
and by Cha & Sembach (2000) are reasonably consistent 
with those determined from the UVES spectrum (Table |6]), 
and the 1996 line profile plotted by Cha & Sembach (2000) 
appears to be very similar to the 2003 UVES profile plot- 
ted in Fig. [1] (though the total equivalent width was smaller 
in 2003). The weak high- velocity components found in the 
UVES spectrum of HD 72127B (similar to those seen to- 
ward several other stars in the Vela region; e.g., Wallerstein 
et al. 1980; Danks & Sembach 1995) do not appear to be 
present in the 1996 spectrum, however. The lines from sev- 
eral of the trace neutral species measured in the 2003 UVES 
spectrum of HD 72127B are significantly stronger than pre- 
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viously reported: the equivalent width of the Na I A5895 hne 
(188 mA) is roughly twice the 95 mA value listed by Hobbs 
et al. (1982), while the equivalent width of the Cai A4226 
line (16 mA) is more than three times the upper limit given 
by Wallerstein & Gilroy (1992). 

3.4 Spatial variations 

Thackeray (1974) first called attention to differences in the 
profiles of the strong, broad Can K lines toward HD 72127A 
and HD 72127B; Hobbs et al. (1982) subsequently noted 
differences in the velocity structures seen in both Call and 
Nal. Higher resolution Call spectra obtained by Cha & Sem- 
bach (2000) indicate that while the overall velocity extent 
of the absorption is very similar toward the two stars, the 
detailed column density distribution (as a function of ve- 
locity) is quite different. Wallerstein et al. (1995b) found 
similar equivalent widths for UV lines of Mg ll and S ll (but 
differences in the strengths of lines from P ll and Ge ll) in 
HST/G160M spectra of the two stars. Possible differences 
in detailed component structure (as seen for Call) could 
not be discerned in those moderate-resolution UV spectra, 
however. 

The UVES spectra of HD 72127A and HD 72127B 
shown in Figs. [1] and [2] provide both a more recent snapshot 
of the differences in Nal and Call and indications of simi- 
larities and differences for lines from several other species. 
Toward HD 72127B, the main component(s) near -f-ll to 
-|-13 km s~^ are stronger in both Nal and Call than those 
toward HD 72127A, but the outlying components generally 
are weaker. The overall structure and relative column den- 
sities are more similar for Nal than for Call - opposite the 
general trend for observations of those two species toward 
other binary systems (Meyer 1990; Watson & Meyer 1996). 
The weak Nal A3302 doublet lines are similar in strength 
toward the two stars, but there are modest differences in 
the profiles of the Ti ll A3383 line and very striking differ- 
ences in the lines from previously undetected Ca I, Cr I, and 
Fel (which are much stronger toward HD 72127B; see next 
section). 

3.5 Trace neutral species toward HD 72127B 

The absorption lines from several trace neutral species are 
unusually strong in the 2003 UVES spectrum of HD 72127B. 
For the iV(H) ~ 2.5x10^" cm"^ estimated above, the col- 
umn density of Nal is somewhat higher than usual (Welty 
& Hobbs 2001), and the column densities of Cal, Cri, and 
Fel are remarkably high (Welty et al. 2003) - especially 
as most of the absorption from those trace neutral species 
is found in the (apparently) single strong component near 
-fll km s~^. The column density of Cal, for example, is 
among the highest known for any Galactic sight line; com- 
parable values generally are seen only for sight lines with 
much higher E{B — V) and/or A'^(H). Absorption from in- 
terstellar Cr I had previously been seen only toward ^ Oph 
(Meyer & Roth 1990), with an equivalent width M^(3578) 
= 0.25±0.06 mA for the strongest line at 3578 A in the 
main (B) components near —15 km s~^, which have N{}i) 
- 14x10^° cm-^ Toward HD 72127B, the Cri A3578 line 
has W(3578) = 5.4±0.7 mA, and five other Cri lines are 
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Table 7. Ratios of trace neutral species (predicted and observed). 



Location Mgl/Nal Cal/Nal Cri/Nal Fel/Nal 



Solar (7000 K) 


+1.55 


+0.08 


-2.43 


+0.92 


Solar (3000 K) 


+0.52 


-1.65 


-2.43 


+0.02 


Solar (100 K) 


+0.60 


-1.29 


-2.43 


+0.10 


Warm (7000 K) 


+1.55 


-1.32 


-3.03 


+0.22 


Warm (3000 K) 


+0.52 


-3.05 


-3.03 


-0.68 


Cold (100 K) 


-0.10 


-4.29 


-4.03 


-1.50 


HD 72127A 


+0.46 


-2.88 


< -1.85 


-1.69 


HD 72127B 




-1.90 


-1.46 


-0.79 


23 Ori SLY 


+0.39 


-3.21 


< -2.82 


-2.07 


C Oph B 


-0.05 


-4.39 


-4.05 


-2.62 


e Ori (-1-3) 


+0.81 


-1.77 






HD 90177 (-40) 




-1.37 




+0.36 


HD 94910 (-54) 




-1.14 




+0.54 


Sk-67 5 (-f261) 




> -0.85 




> +1.05 


Sk-68 52 (-1-307) 




-1.51 




< +0.31 


SN1987A (-1-65) 




-1.57 






SN1987A (-1-165) 




> -0.83 






SN1987A (-1-216) 




-1.74 







Entries are log [Af(X)/Af(Y)]; limits are 3(t. Predicted ratios 
for solar, warm cloud, and cold cloud abundances use the solar 
abundances and representative depletions from Table 5 and 
ionization equilibrium calculations as in eqn. 1 (see also Welty 
et al. 1999a). Predicted values for Cal are for rie = 1.0 cm"''. 
Values for HD 90177 and HD 94910 are from Gnacinski & 
Krogulec (2008); values for Sk-67 5 and Sk-68 52 are from 
Welty &; Crowther (in preparation); values for SN1987A are from 
Welty et al. (1999a). 

seen at greater than 3ct significance - all at the same veloc- 
ity and with equivalent widths consistent with the relative 
/-values (Table [T|. 

Such (relatively) strong absorption from Cal and/or Fel 
has recently been noted in a few other sight lines. Table [7] 
compares the ratios of several trace neutral species Xl/Nal 
- as observed toward HD 72127A, HD 72127B, and several 
other targets - with the values that would be predicted as- 
suming ionization equilibrium under various conditions. The 
last six sight lines in the table exhibit components with even 
more extreme values of Cal/Nal and/or Fel/Nal; the well- 
studied sight lines toward 23 Ori (Welty et al. 1999b) and 
C Oph (e.g., Welty & Hobbs 2001; Welty et al. 2003) are 
included for comparison. The clouds toward HD 90177 and 
HD 94910 are two of the so-called 'CaFe' clouds identified by 
Bondar et al. (2007) and discussed by Gnacinski & Krogulec 
(2008). StiU higher values of Cal/Nal and/or Fel/Nal are 
seen in components toward three targets in the Large Mag- 
ellanic Cloud (Welty et al. 1999a; Wehy & Crowther, in 
preparation). In two of those LMC components, no absorp- 
tion is detected in the Nal D lines, with upper limits on 
iV(Nal) of order 1-4x10^° cm'^ 

Observations of quasar absorption-line systems have be- 
gun to yield detections of weak lines from various trace neu- 
tral species (besides the commonly observed Mg l) in higher 
redshift systems. For example. Si I, Fel, and several other 
neutral species have been detected (along with a modest 
amount of H2) in a sub-Damped Lyman a system &t z ^ 
1.15 (Quast, Reimers, & Baade 2008). Unexpectedly strong 
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absorption from Si I, Cal, and Fel has also been observed 
in a weak Mgll absorber at z ~ 0.45 (D'Odorico 20070. 
Absorption from Nal has not yet been observed in those 
systems, however. 

Consideration of the physical conditions and ionization 
equilibrium in the main cloud toward HD 72127B might 
provide some insight into the high observed abundances of 
the trace neutral species in those various contexts. Waller- 
stein et al. (1995b) found a relatively high fraction of the 
Cl toward HD 72127B to be in the excited fine-structure 
states - indicative of a relatively high thermal pressure 
(tihT) in the gas. High densities would be expected to fa- 
vor the trace neutral species; Cal (and several other neutral 
species) can be enhanced via dielectronic recombination in 
warmer gas. While the width of the Cal component at +13 
km s"^ toward HD 72127A (fe = 0.6 km s"^; Welty et al. 
2003) places an upper limit on the temperature of T < 850 
K for that component (so that dielectronic recombination 
would not be significant there), the main Cal component 
toward HD 72127B may be somewhat broader. Higher res- 
olution spectra of HD 72127B would be needed to place a 
firmer limit on the temperature (and the role of dielectronic 
recombination) in the main component in that sight line. 
The Cal/Call ratio (0.0055) in the main component toward 
HD 72127B, however, is very similar to the ratios seen in 
other Galactic sight lines (e.g., Welty et al. 2003) and is 
only a factor of two higher than the ratio (0.0026) found 
in the main component toward HD 72127A - which would 
seem to suggest that Ca I is not greatly enhanced relative to 
Ca II and that the electron density ne is not especially high. 

Fig. |4]shows the ratios Xl/Nal (for X — Mg, Ca, and 
Fe) that would be predicted for equilibrium between ioniza- 
tion (radiative plus collisional) and recombination (radia- 
tive plus dielectronic), as functions of temperature (as in 
Fig. 11 of Welty et al. 1999a, for example). Those earlier 
calculations (which assumed 'warm cloud' depletions) have 
been adjusted here for total gas phase elemental abundances 
characterized also by the patterns seen in the solar system 
(Lodders 2003) and in interstellar clouds exhibiting 'cold 
cloud' depletions (e.g., Savage & Sembach 1996; Welty et 
al. 1999b; see Table [5]). The first six lines of Table[7]hst the 
predicted ratios for several representative combinations of 
depletion and temperature. For Mg, Cr, and Fe - for which 
the first ions X ll should be dominant - the predicted ratios 
are given by 



Xl 
Nal 



Ax (5x ax (aNaWe + FNa + CNa^e) 



(1) 



ANa ^Na ONa (axrie + Fx + CxMe) 

where Ax and (5x are the solar abundance and depletion 
of element X, respectively. The calculation for calcium is 
slightly more complicated, as Cam must also be consid- 
ered. The adopted photoionization rates (Fx) are taken in 
most cases from Pequignot & Aldrovandi (1986), assuming 
the WJl radiation field (de Boer et al. 1973); the collisional 

Based on comparisons with Galactic sight lines, D'Odorico 
(2007) infers log[7V(H)] ~ 20.5-21.5 - and thus a very low metal- 
licity - for that absorber. The Galactic sight lines, however, are 
typically characterized by significant depletions of the various re- 
fractory elements. The Af(H) in the weak Mgll system (where the 
depletions are likely much less severe) may well be much lower, 
and the metallicity correspondingly higher. 



a. 



4 
3 
2 
1 

-1 
-2 



Mg I/Na I 




log(T) 



^ -2 
Id 

Z. -3 



log(T) 



-3 



-4 



3 





03 

o 



log(T) 

Figure 4. Predicted ratios Xl/Nal - for X = Mg, Ca, and Fe - 
as functions of temperature. The calculations assume equilibrium 
between ionization (radiative and collisional) and recombination 
(radiative and dielectronic), warm cloud depletions, and = 
0.001, 0.01, 0.1, or 1.0 cm~^ (lower to upper curves in each panel). 
The lower partial curves (from 100 to 300 K) show the ratios for 
cold cloud depletions. The increases at higher temperatures are 
due to dielectronic recombination to Mg I, Ca I, and Fe I. The lower 
set of curves for Fe I shows the predictions for the photoionization 
and recombination rates calculated by Nahar et al. (1997). The 
letters 'S' and 'B' mark the ratios observed for the 23 Ori SLV 
components and the C, Oph B components, respectively; the bar at 
right shows the range for each ratio observed in the local Galactic 
ISM. 
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ionization rates (cx) are taken from Shull & Van Steenberg 
(1982). The total (radiative plus dielectronic) recombination 
rate coefficients (ax) are calculated from the parameters 
given by Aldrovandi & Pequignot (1973, 1974), Shull & Van 
Steenberg (1982), and/or Pequignot & Aldrovandi (1986). 
The lower set of curves for Pel shows the ratios predicted 
using the photoionization and (total) recombination rates 
computed by Nahar, Bautista, & Pradhan (1997). For Cri, 
Pcr = 8 X 10"^° s^^ was estimated by Meyer & Roth (1990), 
and acr has been set to 6 x 10~^^ cm^s"^ (for T = 100 K), 
similar to the values determined for radiative recombina- 
tion to a number of other trace neutral species (Pequignot 
& Aldrovandi 1986); the temperature dependence of acr is 
assumed to be the same as for ctNa (i.e., with no significant 
contribution from dielectronic recombination for T < 12000 
K; see comments below, however). If the recombination is 
dominated instead by charge exchange with large molecules 
or small grains (Lepp et al. 1988; Weingartner & Draine 
2001; Liszt 2003; Welty et al. 2003), and if the resuhing 
neutrals do not stick to the grains (sticking parameter s = 
0), then the predicted ratios change by at most 0.1 dex. 

The predicted ratios for Mgl, Cal, and Fel shown in 
Fig. |4] exhibit some common trends. At temperatures below 
about 3000 K, the ratios reflect primarily photoionization 
and radiative recombination - and are not very sensitive 
to the overall strength of the radiation field, the electron 
density rie, or the temperature. Because Ca ll is often a trace 
species, however, the Cal/Nal ratio does depend on rie. At 
somewhat higher temperatures, all the ratios increase due 
to dielectronic recombination (which does not significantly 
affect Nal until much higher temperatures) - starting at 
about 4000 K, 3000 K, and 5000 K for Mgi, Cai, and Fei, 
respectively (Aldrovandi & Pequignot 1973, 1974; Shull & 
Van Steenberg 1982). 

For the main cloud(s) toward Oph, the observed ra- 
tios for Mgl, Cal, and Cri are in good agreement with the 
predicted cold cloud values (consistent with the depletions 
derived from the corresponding dominant species and the 
temperature inferred from H2 rotational excitation). The 
ratio for Fel is lower than the predicted cold cloud value 
by more than a factor of 10, however [but is not as deficient 
relative to the value predicted using the Nahar et al. (1997) 
rates]. Meyer & Roth (1990) remarked that the electron den- 
sity estimated toward Oph from the Cri/Crll ratio was 
similar to the values obtained from other neutral/first ion 
ratios; Welty et al. (2003) noted that it is not uncommon for 
Fe I to be somewhat weaker than expected, relative to other 
trace neutral species. For the 'strong low-velocity' clouds to- 
ward 23 Ori (where some trace neutral species are relatively 
strong and the depletions are intermediate between the rep- 
resentative warm and cold cloud values; Welty et al. 1999b) 
and for the main cloud(s) toward HD 72127A (with warm 
cloud depletions), the observed ratios for Mgl and Cal are 
closer to the predicted warm cloud values for T ~ 3000 K 
(although the gas is likely much cooler than that in both 
cases), but the ratios for Fel again are low by factors of 10- 
25. For the main cloud(s) toward HD 72127B, the observed 
ratios for both Cal and Fel are about a factor of 10 higher 
than those toward HD 72127A. The ratio for Cal thus is 
closer to the value predicted for solar relative abundances 
than to that predicted for warm cloud abundances (for T 
^ 3000 K), while the ratio for Fel is consistent with the 



predicted warm cloud value at that T (i.e., low compared to 
Cal); the observed ratio for Cri, however, is a factor of 10 
higher than the predicted solar value. 

As calcium, chromium, and iron usually are much more 
severely depleted into dust than is sodium, an enhancement 
of the neutral species Cal, Cri, and Fel, relative to Nal, 
could be due to much less severe depletions. As noted above, 
dielectronic recombination can increase the abundances of 
Mgl, Cal, and Fel when the temperature exceeds about 
4000 K, 3000 K, and 5000 K, respectively. The cloud at -1-65 
km s"^ toward SN1987A is relatively cold [as the Nal D 
lines at that velocity exhibit resolved hyperfine structure, 
with b ~ 0.3 km s"^ (Pettini & Gillingham 1988; Welty & 
Crowther, in preparation)], but the ratios of various dom- 
inant species suggest that there is essentially no depletion 
there (Welty et al. 1999a). The observed Cal/Nal ratio is 
within a factor of 2 of the predicted value for solar rel- 
ative abundances at 100 K (i.e., a cold, dust-free cloud), 
though it is also consistent with the values predicted for 
somewhat warmer gas. For the component at -1-3 km s~^ to- 
ward e Ori, the 6- values for Na I, Ca I, and Ca 11 - all 0.36-0.38 
km s~^ - indicate that that cloud also is fairly cold (Welty et 
al. 1994, 1996, 2003). The slightly elevated Mgl/Nal ratio 
thus may reflect less severe depletion of magnesium there 
(and not dielectronic recombination). The still higher val- 
ues of Cal/Nal and/or Fel/Nal seen in several other cases, 
however, would seem to imply some combination of mild 
(or negligible) depletions and dielectronic recombination in 
warmer gas (e.g., Gnacihski & Krogulec 2008 j3- though the 
Cal/Call ratios generally are not unusually high. 

The strong observed absorption from Cal, Cri, and 
Fel toward HD 72127B - relative to Nal and to the much 
weaker absorption from those species toward HD 72127A 
- also may suggest some combination of mild depletions 
and higher temperature. The similarities in the strengths 
of the 5780 A DIB, the E{B - V) colour excesses, and the 
equivalent widths of the weak Mgll A1239,1240 lines toward 
HD 72127A and HD 72127B suggest that the total hydrogen 
column densities are similar for the two sight lines. The rel- 
atively small differences in the column densities of Na I and 
Ti II in the two sight lines (both overall and in the main com- 
ponents near -1-12 km s~^) and in the Cal/Call ratios (in the 
main components) then suggest that differences in ionization 
and depletions could account for only part of the enhance- 
ment of Cal, Cri, and Fel. The higher values for Cal/Nal 
and Fel/Nal toward HD 72127B could be explained, how- 
ever, if the Cal and Fel were enhanced via dielectronic re- 
combination, in gas with warm cloud depletions and T > 



^ Gnacinski & Krogulec (2008) assume that there is no depletion 
in the CaFe clouds - taking non-detections of CH to imply a lack 
of both H2 and dust. We note, however, that there are many sight 
lines where CH is not detected (to similar limits), but which none 
the less have E{B — V) greater than zero, detections of H2 (at 
column densities below the regime of self-shielding), and evidence 
from dominant species for at least moderate depletions. The high 
temperatures ( > 8000 K) derived for the CaFe clouds presumably 
reflect the ability of dielectronic recombination to produce the 
observed abundances of Cal and Fe I in such warm gas. Allowance 
for depletions and for colder temperatures, however, could reduce 
the substantial discrepancies between their models and the Nal 
and Kl observed for the 'non-CaPe' clouds. 
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5000 K. We conjecture that dielectronic recombination to 
Cri might then account for the high Cri/Nal ratio toward 
HD 72127B - and we would predict some enhancement of 
Mgl (which is not unusually strong toward HD 72127A), 
but no enhancement of Si I (for which dielectronic recom- 
bination becomes important at somewhat higher tempera- 
tures If dielectronic recombination is responsible for the 
strong absorption from those trace neutral species toward 
HD 72127B, then there would be a significant difference in 
temperature for the main components toward HD 72127A 
and HD 72127B - over a distance of less than 2200 AU. 
We note, however, that the Cal/Call ratio often seems to 
imply a higher rie than other such trace/dominant ratios - 
even when the Cal 6- values place strong constraints on T 
(Welty et al. 2003) - so that high Cal abundances are not 
necessarily due to dielectronic recombination in warm gas. 

3.6 Location of gas 

The location of the main absorption components a.t v ^ 
-1-10 to -1-13 km s~^ (ulsr ~ —3 to km s"'^) toward 
HD 72127AB - with depletions apparently similar to those 
found in typical warm, diffuse disc clouds - is not known. 
The outlying components, with higher LSR velocities and 
much milder depletions, are very likely associated with the 
Vela SNR, at a distance of about 250±30 pc (Cha et al. 
1999). Some absorption at low LSR velocities, however, is 
seen for stars closer than 200 pc - and appears to be located 
at distances less than about 130 pc (Cha et al. 2000). The 
foreground absorption is 'patchy and inhomogeneous' (Cha 
et al. 2000), and the highest column density of Nal seen 
toward any star foreground to the SNR is about 7.2x10^^ 
cm~^ (the V = 12.3 km s~^ component toward HD 72232, 
at less than 115 pc) - a factor of 3-5 lower than the values 
obtained for the main components toward HD 72127AB. 
The relatively weak Call absorption toward HD 72232 (15 
mA; Cha & Sembach 2000) implies a Can column density 
of about 2x10^^ cm~^; the corresponding Call/Nal ratio, 
about 0.3, is similar to the value found for the main compo- 
nent toward HD 72127A (Table[21). The much higher column 
densities and the continuing variations in both column den- 
sity and velocity suggest that the bulk of the material in the 
main components toward HD 72127AB is associated with 
the SNR, but there is very likely some contribution from 
foreground material (at the same velocity) as well. 



4 SUMMARY / CONCLUSIONS 

New, moderately high resolution (FWHM = 4.5-4.9 
km s~^) optical spectra of Nal and Call absorption toward 
HD 72127AB have provided additional evidence for both 
spatial and temporal variations in the complex interstellar 
absorption along the two sight lines. If the spatial differences 
are associated with material in the Vela SNR, they occur on 

* After the initial submission of this paper, we received prelim- 
inary, partial calculations of the total recombination rate coef- 
ficients to Cri (S. Nahar, private communication), which both 
confirm our assumptions regarding the recombination at T < 
1000 K and suggest that dielectronic recombination may indeed 
be significant below 10"* K for Cr I. 



scales of about 1100 AU. Fits to the absorption-line profiles 
seen in high-resolution (FWHM ~ 3.5 km s~^) UV spectra 
of HD 72127A obtained with the HST GHRS have yielded 
abundances for a number of species (e.g.. Sill, Fell, Nill, 
and Znll) in the various velocity components discernible 
in the spectra. The main components at low LSR veloci- 
ties have depletions similar to those found in warm, diffuse 
disc clouds; the generally weaker components at higher ve- 
locities have much milder depletions, more similar to those 
found in halo clouds. Similarities in E{B — V) and in the 
equivalent widths of the 5780 A diffuse interstellar band and 
the UV lines of Mg ll and S ll suggest that the total hydro- 
gen column densities in the two sight lines are similar, with 
log[A'^(H) (cm~^)] ~ 20.40. Similarities in the column den- 
sities of Ti II and Na I suggest that the depletions and ion- 
ization in the main components also are comparable (within 
factors of about 2). Several other trace neutral species - 
Cal, Cri, and Fel - are much stronger toward HD 72127B, 
however. In particular, the column density of Cri is about 
30 times the value found for the main components toward 
C, Oph (the only other sight line in which Cri has been de- 
tected) - even though A'^(H) is a factor of about 6 lower 
toward HD 72127B. While an earlier high-resolution spec- 
trum of Cal toward HD 72127A indicates that the main 
component in that sight line has T < 900 K, the strong 
absorption from Cal and Fel toward HD 72127B suggests 
that the main component there is much warmer - T > 5000 
K - with dielectronic recombination largely responsible for 
the enhanced abundances of those neutral species. We con- 
jecture that dielectronic recombination may also be respon- 
sible for the enhanced Cri absorption toward HD 72127B, 
and predict that strong absorption will be found for Mg I as 
well. 

It would be very useful to obtain high-resolution 
UV spectra of HD 72127B: Sill, Cm, Feii, and Znii 
[to determine Af(H) and depletions and to compare with 
HD 72127A]; Mgl and Si I (to obtain constraints on the 
effect of dielectronic recombination) ; and C I (also toward 
HD 72127A; to estimate thermal pressures and densities in 
the various components). Higher resolution optical spectra 
of the trace neutral species would provide more stringent 
constraints on the temperature in the main component. De- 
termination of the recombination behaviour of Cm, as a 
function of temperature, would aid in understanding the 
high abundance of Cr I and the role of dielectronic recombi- 
nation. Continued monitoring of the interstellar absorption 
toward these two stars could provide insights into the struc- 
ture and properties of interface regions between warm and 
cool gas. 
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